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WO 99/50773 PCT/US99/05463 

Portfolio Optimization by Means of Resampled Efficient Frontiers 

Technical FipJrf 

The present invention relates to a method for selecting a portfolio of tangible or 
5 intangible assets subject to optimization criteria yielding an enhanced mean-variance 
efficiency. 

Background of the Invention 
Managers of assets, such as portfolios of stocks, projects in a firm, or other 

10 assets, typically seek to maximize the expected or average return on an overall 

investment of funds for a given level of risk as defined in terms of variance of return, 
either historically or as adjusted using techniques known to persons skilled in portfolio 
management. Alternatively, investment goals may be directed toward residual return 
with respect to a benchmark as a function of residual return variance. Consequently, the 

15 terms "return" and ''variance," as used in this description and in any appended claims, 
may encompass, equally, the residual components as understood in the art. The capital 
asset pricing model of Sharpe and Lintner and the arbitrage pricing theory of Ross are 
examples of asset evaluation theories used in computing residual returns in the field of 
equity pricing. Alternatively, the goal of a portfolio management strategy may be cast as 

20 the minimization of risk for a given level of expected return. 

The risk assigned to a portfolio is typically expressed in terms of its variance Op 
stated in terms of the weighted variances of the individual assets, as: 



25 

where w A is the relative weight of the i-th asset within the portfolio, 

30 is the co variance of the i-th and j-th assets, is their correlation, and o { is the standard 
deviation of the i-th asset The portfolio standard deviation is the square root of the 
variance of the portfolio. 
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Following the classical paradigm due to Markowitz. a portfolio may be 
optimized, with the goal of deriving the peak average return for a given level of risk and 
any specified set of constraints, in order to derive a so-called "mean-variance (MV) 
efficient* portfolio using known techniques of linear or quadratic programming as 
appropriate. Techniques for incorporating multiperiod investment horizons are also 
known in the an. As shown in FIG. 1 A. the expected return p for a portfolio may be 
plotted versus the portfolio standard deviation a, with the locus of MV efficient 
portfolios as a function of portfolio standard deviation referred to as the "MV efficient 
frontier/* and designated by the numeral 10. Mathematical algorithms for deriving the 
MV efficient frontier are known in the art. 

Referring to FIG. 1B ? a variation of classical Markowitz MV efficiency often 
used is benchmark optimization. In this case, the expected residual return a relative to a 
specified benchmark is considered as a function of residual return variance u>. defined as 
was the portfolio standard deviation a but with respect to a residual risk. An investor 
with portfolio A desires to optimize expected residual return at the same level u> A of 
residual risk. As before, an efficient frontier 10 is defined as the locus of all portfolios 
having a maximum expected residual return a of each of all possible levels of portfolio 
residual risk. 

Known deficiencies of MV optimization as a practical tool for investment 
management include the instability and ambiguity of solutions. It is known that MV 
optimization may give rise to solutions which are both unstable with respect to small 
changes (within the uncertainties of the input parameters) and often non-intuitive and 
thus of little investment sense or value for investment purposes and with poor out-of- 
sample average performance. These deficiencies are known to arise due to the propensity 
of MV optimization as "estimation-error maximizers,'* as discussed in R. Michaud. "The 
Markowitz Optimization Enigma: Is Optimized Optimal?" Financial Analysts Journal 
(1989), which is herein incorporated by reference. In particular. MV optimization tends 
to overweight those assets having large statistical estimation errors associated with large 
estimated returns, small variances, and negative correlations, often resulting in poor ex- 
post performance. 
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Summary of the Invention 
In accordance with one aspect of the invention, in one of its embodiments, there 
is provided a method for evaluating an existing or putative portfolio having a plurality of 
assets. The existing portfolio is of the kind having a total portfolio value, where each 
asset has a value forming a fraction of the total portfolio value, each asset has a defined 
expected return and a defined standard deviation of return, and each asset has a 
covariance with respect to each of every other asset of the plurality of assets. The method 
has the steps of: 

a. computing a mean-variance efficient frontier based at least on input data 
characterizing the defined expected return and the defined standard deviation of return of 
each of the plurality of assets: 

b. indexing a set of portfolios located on the mean variance efficient frontier; 

c. resampling a plurality of simulations of input data statistically consistent with 
the defined expected return and the defined standard deviation of return of each of the 
plurality of assets: 

d. computing a simulated mean-variance efficient portfolio for each of the 
plurality of simulations of input data; 

e. associating each simulated mean-variance efficient portfolio with a specified 
portfolio of the set of indexed portfolios for creating a set of identical-index-associated 
mean-variance efficient portfolios; 

f. establishing a statistical mean for each set of identical- index-associated mean- 
variance efficient portfolios, the plurality of statistical means defining the resampled 
efficient frontier: and 

g. comparing the existing portfolio with a portfolio from the resampled efficient 
frontier characterized by a specified risk objective. 

In accordance with an alternate embodiment of the present invention, an index- 
associated resampled fuzzy region may be identified in the risk/return plane 
corresponding to portfolios having a norm, with respect to a corresponding portfolio on 
the resampled efficient frontier, that is less than or equal to a value corresponding to a 
specified confidence level. A current portfolio may be assessed with respect to required 
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optimization based at least on the value of the norm of the current portfolio relative to 
the index-associated resampled fuzzy region. 

Brief Description of the Drawings 
5 The invention will be more readily understood by reference to the following description, 
taken with the accompanying drawings, in which: 

FIGS. 1A depicts the prior an principle of calculating an efficient frontier of 
maximum expected return for each given level of portfolio risk: 

FIGS. IB depicts the prior art principle of calculating an efficient frontier of 
10 maximum expected residual return for each given level of portfolio residual risk: 

FIG. 2 displays a set of statistically equivalent portfolios within the risk/return 

plane: 

FIG. 3 displays statistically equivalent portfolios within the risk/return plane 
corresponding to three particular risk rankings on the efficient frontier: minimum 
15 variance, maximum return, and a middle return portfolio; 

FIG. 4 displays a region within the risk/return plane corresponding to statistically 
equivalent portfolios for a given level of probability; 

FIG. 5 shows the resampled efficient frontier plotted in the risk/return plane in 
accordance with a preferred embodiment of the present invention: and 
20 FIG. 6 displays a fuzzy region of statistically equivalent portfolios in the 

risk/return plane relative to resampled efficiency for two particular risk rankings for a 
given level of probability: minimum variance and middle return. 

Detailed Description of Preferred Embodiments 
25 The foregoing features of the invention will be more readily understood by 

reference to the following detailed description taken with the accompanying drawings. 

In accordance with various embodiments of the present invention, it is 
recognized that MV optimization is a statistical procedure, based on estimated returns 
subject to a statistical variance, and that, consequently, the MV efficient frontier, as 
30 defined above, is itself characterized by a variance. For a given set of input data 

characterizing the universe of assets to be considered to comprise a portfolio, and for a 
given set of constraints on the makeup of the portfolio, an MV efficient frontier may be 

4 
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calculated, using standard techniques discussed above. Since the input data are of a 
statistical nature (i.e., characterized by means with associated variances and other 
statistical measures), the input data may be resampled. by simulation of optimization 
input parameters in a manner statistically consistent with the first set of data, as 
5 described, for example, by J. Jobson and B. Korkie, "Estimation for Markowuz Efficient 
Portfolios," Journal of Portfolio Mangement . (1981), which is herein incorporated by 
reference. Other resampling methods are within the scope of the invention as described 
herein, and as claimed in arry appended claims. 

Based on the resampled data, the MV efficient portfolios may be recalculated. 

10 subject to the same constraints as applied in the initial solution. In this manner, as shown 
in Fig. 2, a set 12 of statistically equivalent MV efficient portfolios may be calculated. 
By iterating this procedure, a large MV efficient "statistical equivalence'" set of 
portfolios, in the expected return - portfolio variance space, may be generated. In 
particular, an existing portfolio may be found to be consistent, in a statistical sense, with 

15 MV efficiency and thus not to require optimization, thereby potentially and 

advantageously saving transaction costs associated with revision of a portfolio. Referring 
to Fig. 3, statistically equivalent portfolios within the risk/return plane are shown 
corresponding to three particular risk rankings on the efficient frontier: namely, 
minimum variance 14. maximum return 16. and a middle return portfolios 18. 

20 One method of establishing statistical equivalence is described in reference to 

FIG. 4. In accordance with one sample acceptance procedure, a sample acceptance area 
20 is defined under efficient frontier 10. The 100( l-a)% sample acceptance region 20 
includes, on average, the fraction 1-oc of Monte-Carlo estimated resampled portfolios, 
where a is the probability of a Type I error, i.e., the probability of rejecting the null 

25 hypothesis when it is true, which are statistically equivalent to a portfolio denoted as 

CURRENT. The loci corresponding, in the risk/return plane, to particular assets included 
within the universe from which the sample portfolio is drawn are indicated in the figure. 
Thus, for cc=0.1, the 90% sample acceptance region 20 is obtained. Sample acceptance 
region 20 may be obtained, in practice, by a procedure in which the area under efficient 

30 frontier 10 is divided into mutually exclusive column (or row) rectangles that include all 
the simulated portfolios. The base of the rectangle is defined as the minimum return 
point that contains the fraction 1-oc of the simulated portfolios within the rectangle. The 
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midpoints of the rectangles are connected to obtain curve 22 which constitutes the lower 
boundary of sample acceptance area 20. Irregularities in the shape of curve 22 are 
attributable in part to the small number of assets included in the example. 

A disadvantage of the foregoing method is that all portfolios which are 
5 statistically equivalent (in terms of the specified criterion) appear within sample 

acceptance region 20 without accounting for desired risk objectives or other constraints. 
Referring now to Fig. 3, a method is described whereby, in accordance with a preferred 
embodiment of the present invention, all resampled portfolios within the risk/return 
plane may be associated, many-to-one, with particular portfolios on MV efficient frontier 

10 10. Various criteria may be applied in associating portfolios with those on the MV 

efficient frontier, and all such associations are within the scope of the present invention. 
As one example, each efficient frontier portfolio (i.e.. each point on efficient frontier 10) 
may be identified by its relative return rank. Similarly, the efficient frontier portfolios 
may be ranked by their variance, the maximum variance corresponding to the maximum 

15 return, the rankings by risk or return similarly mapping onto one another uniquely. Thus, 
for example, the minimum variance portfolio 14 might have the lowest rank relative to 
the other efficient portfolios of efficiency frontier 10. Similarly, maximum average 
return portfolio 16 has the highest average return rank in each simulated efficient 
frontier. Similarly, any other simulated portfolio is rank associated with a particular 

20 efficient frontier portfolio. The sparsely clustered portfolios 18 shown in the figure 
correspond to the middle' ranked efficient portfolio. In practice, the shape of the rank- 
associated regions varies in dependence upon the position of the portfolio on the MV 
efficient frontier. 

It is not necessary, however, that the association with efficient frontier portfolios 
25 be by rank, and particular portfolios on the MV efficient frontier may be indexed, and 
thus index -associated, each with a set of statistically equivalent efficient portfolios lying 
below the efficient frontier. 

In accordance with an alternate embodiment of the invention, indexing of the set 
of MV efficient portfolios for associating the resampled portfolios is to associate with 
30 each MV efficient portfolio a "lambda value," defining the risk/return preference, with 
respect to which the quantity <|)=a 2 -Aji is minimized, where a 2 is the variance of each 
portfolio and p is the expected return of each portfolio of the set of portfolios located on 
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the mean variance efficient frontier. The parameter X assumes a value between zero and 
infinity. 

Once the ensemble of index-associated MV efficient portfolios has been 
associated, usual statistical measures of the ensemble may be derived. These measures 
include, without limitation, the averages, standard errors, and t-statistics of the average 
of the portfolio weights of the rank-associated simulated efficient portfolios. Referring 
now to Fig. 5, an average of index-associated MV efficient portfolios may be defined, in 
accordance with preferred embodiments of the present invention, the average of index- 
associated MV efficient portfolios being referred to as a "resampled-efficient portfolio." 
The average may be determined with respect to any of a variety of parameters, and, in 
accordance with a preferred embodiment, it is with respect to the vector average of the 
associated portfolios. The locus 40 of resampled-efficient portfolios is referred to as the 
"resampled efficient frontier." The resampled-efficient portfolio and its associated 
statistics may be applied as a statistical measure for portfolio analysis, as further 
described below. Its application, as a choice for portfolio selection, advantageously 
removes, by definition, the ^outlier" portfolios which strongly depend on values of a 
particular set of inputs and improves out-of-sample performance, on average. Statistical 
procedures for portfolio analysis and revision, and performance benefits based on these 
concepts, in accordance with embodiments of the present invention, are also described in 
detail in the book, R. Michaud, Efficient Asset Management, (Harvard Business School 
Press. 1998), which is hereby incorporated by reference, as are all references cited in the 
text of the book. 

In accordance with other embodiments of the invention, the portfolio 
optimization problem may be cast in terms of fuzzy sets and the application of fuzzy 
logic. Referring now to Fig. 6. the "similarity" of a particular resampled efficient frontier 
portfolio (whether risk-indexed or otherwise indexed) to a corresponding resampled 
efficient portfolio may be defined by a "norm/" with the norm having the usual 
properties of a distance function as known to persons skilled in the mathematical arts. 
Risk-indexing will be assumed in the following description solely for purposes of 
coherence of presentation. In a preferred embodiment, the norm may be a Euclidean 
vector distance, or U norm. The distribution of the values of the norms of risk-index 
associated resampled portfolios relative to the associated resampled efficient frontier 
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portfolio defines the distance criterion and fuzzy or similarity resampled efficient region 
(again, relative to the associated resampled efficient frontier portfolio) at a given level of 
probability. A fuzzy region, designated generally by numeral 60. corresponds to the 
region in the risk/return plane at probability level 0.9 using the Euclidean norm for the 
5 minimum variance and middle resampled efficient risk-indexed portfolios, each taken 
relative to the resampled efficient frontier. In particular, the distance criterion for any 
portfolio P is taken to be the relative variance for portfolio P, 

(p-p 0 y * s * (p-p q ) . 

where P-P 0 is the difference vector of portfolio weights with respect to P nf the 
corresponding index- associated portfolio on the resampled efficient frontier, and. S is the 

10 input return covariance matrix (with the superscript Y denoting the transpose of the 
difference vector). The norm is taken in the space of portfolio vectors (i.e., "portfolio 
space"). The distance criterion defined for probability level 1-a is the value of the norm 
of the risk-indexed resampled portfolios corresponding to the 100(1 -a) percentile value 
of all risk-indexed resampled portfolios sorted from small to large. An example of the 

15 use of a confidence region resampling test statistic relative to the resampled efficient 
frontier is provided at page 73 of Michaud, Efficient Asset Management. Portfolios 
within a risk-indexed resampled fuzzy region at a given level of probability may be said 
to be sufficiently similar to the risk-indexed resampled efficient portfolio to be 
considered resampled efficient and not to require optimization. Particular probability 

20 levels and Lp norms may be specified as a matter of design choice, within the scope of 
the present invention. 

In an alternative embodiment, the disclosed method for evaluating an existing or 
putative portfolio may be implemented as a computer program product for use with a 
computer system. Such implementation may include a series of computer instructions 

25 fixed either on a tangible medium, such as a computer readable medium {e.g., a diskette, 
CD-ROM, ROM ? or fixed disk) or transmittable to a computer system, via a modem or 
other interface device, such as a communications adapter connected to a network over a 
medium. The medium may be either a tangible medium (e.g., optical or analog 
communications lines) or a medium implemented with wireless techniques {e.g., 

30 microwave, infrared or other transmission techniques). The series of computer 
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instructions embodies all or part of the functionality previously described herein with 
respect to the system. Those skilled in the art should appreciate that such computer 
instructions can be written in a number of programming languages for use with many 
computer architectures or operating systems. Furthermore, such instructions may be 
5 stored in any memory device, such as semiconductor, magnetic, optical or other memory 
devices, and may be transmitted using any communications technology, such as optical, 
infrared, microwave, or other transmission technologies. It is expected that such a 
computer program product may be distributed as a removable medium with 
accompanying printed or electronic documentation (e.g., shrink wrapped software). 

10 preloaded with a computer system (e.g., on system ROM or fixed disk), or distributed 
from a server or electronic bulletin board over the network (e.g., the Internet or World 
Wide Web). Of course, some embodiments of the invention may be implemented as a 
combination of both software (e.g., a computer program product) and hardware. Still 
other embodiments of the invention are implemented as entirely hardware, or entirely 

1 5 software (e.g., a computer program product). 

The described embodiments of the invention are intended to be merely exemplary 
and numerous variations and modifications will be apparent to those skilled in the art. 
All such variations and modifications are intended to be within the scope of the present 
invention as defined in the appended claims. 
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WE CLAIM: 

1. A method tor selecting a value of portfolio weight for each of a specified 
plurality of assets of an optimal portfolio, the value of portfolio weight chosen from 
values between zero and unity, each asset having a defined expected return and a defined 
standard deviation of return, each asset having a covariance with respect to each of every 
other asset of the plurality of assets, the method comprising: 

a. computing a mean-variance efficient frontier based at least on input data 
characterizing the defined expected return and the defined standard deviation of return of 
each of the plurality of assets; 

b. indexing a set of portfolios located on the mean variance efficient frontier; 

c. resampling a plurality of simulations of input data statistically consistent with 
the defined expected return and the defined standard deviation of return of each of the 
plurality of assets: 

d. computing a mean-variance efficient portfolio for each of the plurality of 
simulations of input data; 

e. associating each simulated mean-variance efficient portfolio with a specified 
portfolio of the set of indexed portfolios for creating a set of identical-index-assoctated 
mean-variance efficient portfolios; 

f. establishing a statistical mean for each set of identical-index-associated mean- 
variance efficient portfolios, the plurality of statistical means defining the resampled 
efficient frontier; and 

g. selecting a portfolio weight for each asset from the resampled efficient 
frontier according to a specified risk objective. 

2. A method according to claim 1, wherein the step of indexing the set of portfolios 
includes associating a rank with each portfolio of the set of portfolios located on the 
mean variance efficient frontier. 

3. A method according to claim 1, wherein the step of indexing the set of portfolios 
includes associating a lambda value for which 4>=a 2 -A.|u is minimized, where o z is the 
variance of each portfolio and jli is the expected return of each portfolio of the set of 
portfolios located on the mean variance efficient frontier. 

4. A method for evaluating an existing portfolio having a plurality of assets, the 
existing portfolio having a total portfolio value, each asset having a value forming a 

!() 
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fraction of the total portfolio value, each asset having a defined expected return and a 
defined standard deviation of return, each asset having a covariance with respect to each 
of every other asset of the plurality of assets, the method comprising: 

a. computing a mean-variance efficient frontier based at least on input data 

5 characterizing the defined expected return and the defined standard deviation of return of 
each of the plurality of assets; 

b. associating a rank with each portfolio of assets associated with the mean 
variance efficient frontier in such a manner as to create a plurality of ranked portfolios: 

c. resampling a plurality of simulations of input data statistically consistent with 
10 the defined expected return and the defined standard deviation of return of each of the 

plurality of assets; 

d. computing a mean- variance efficient portfolio for each of the plurality of 
simulations of input data: 

e. associating each simulated mean-variance efficient portfolio with a specified 
5 portfolio of the set of indexed portfolios for creating a set of identical-index-associated 

mean- variance efficient portfolios; 

f. establishing a statistical mean for each set of identical-index-associated mean- 
variance efficient portfolios, the plurality of statistical means defining the resampled 
efficient frontier: and 

0 g. comparing the existing portfolio with a portfolio from the resampled efficient 

frontier characterized by a specified risk objective. 

5. A method for evaluating an existing portfolio having a plurality of assets, the 
existing portfolio having a total portfolio value, each asset having a value forming a 
fraction of the total portfolio value, each asset having a defined expected return and a 
5 defined standard deviation of return, each asset having a covariance with respect to each 
of every other asset of the plurality of assets, the method comprising: 

a. computing a mean-variance efficient frontier based at least on input data 
characterizing the defined expected return and the defined standard deviation of return of 
each of the plurality of assets: 

b. associating a rank with each portfolio of assets associated with the mean 
variance efficient frontier in such a manner as to create a plurality of ranked portfolios: 

11 
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c. resampling a plurality of simulations of input data statistically consistent with 
the defined expected return and the defined standard deviation of return of each of the 
plurality of assets; 

d. computing a mean-variance efficient portfolio for each of the plurality of 
5 simulations of input data: 

e. associating each simulated mean-variance efficient portfolio with a specified 
portfolio of the set of indexed portfolios for creating a set of identical-index-associated 
mean- variance efficient portfolios; 

f. establishing a statistical mean for each set of identical-index-associated mean- 
10 variance efficient portfolios, the plurality of statistical means defining the resampled 

efficient frontier; and 

g. associating a norm having a value characterizing a similarity of each portfolio 
with respect to a corresponding portfolio on the resampled efficient frontier; 

h. defining a index-associated resampled fuzzy region corresponding to 

15 portfolios having a norm value less than or equal to a value associated with a specified 
confidence level: and 

i. evaluating whether a current portfolio requires optimization based at least on 
a norm of the current portfolio in comparison with the index-associated resampled fuzzy 
region. 

6. A computer program product for use on a computer system for selecting a value 
of portfolio weight for each of a specified plurality of assets of an optimal portfolio, the 
value of portfolio weight chosen from values between zero and unity, each asset having a 
defined expected return and a defined standard deviation of return, each asset having a 
covariance with respect to each of every other asset of the plurality of assets, the 
computer program product comprising a computer usable medium having computer 
readable program code thereon, the computer readable program code including: 

a. program code for computing a mean-variance efficient frontier based at least 
on input data characterizing the defined expected return and the defined standard 
deviation of return of each of the plurality of assets: 

b. a sequencer for indexing a set of portfolios located on the mean variance 
efficient frontier: 
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c. a routine for resampling a plurality of simulations of input data statistically 
consistent with the defined expected return and the defined standard deviation of return 
of each of the plurality of assets; 

d. program code for computing a mean-variance efficient portfolio for each of 
5 the plurality of simulations of input data: 

e. program code for associating each simulated mean-variance efficient portfolio 
with a specified portfolio of the set of indexed portfolios for creating a set of identical- 
index-associated mean-variance efficient portfolios: 

f. a module for establishing a statistical mean for each set of identical-index- 
10 associated mean- variance efficient portfolios, the plurality of statistical means defining 

the resampled efficient frontier; and 

g. program code for selecting a portfolio weight for each asset from the 
resampled efficient frontier according to a specified risk objective. 
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